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Abstract 

It is shown that the so-called generic instabilities that appear in the framework of rel- 
ativistic linear irreversible thermodynamics, describing the fluctuations of a simple fluid 
close to equilibrium, arise due to the coupling of heat with hydrodynamic acceleration 
which appears in Eckart's formalism of relativistic irreversible thermodynamics. Further, 
we emphasize that such behavior should be interpreted as a contradiction to the postulates 
of linear irreversible thermodynamics (LIT), namely a violation of Onsager's hypothesis 
on the regression of fluctuations, and not as fluid instabilities. Such contradictions can be 
avoided within a relativistic linear framework if a Meixner-like approach to the phenomeno- 
logical equations is employed. 
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I. INTRODUCTION 



Relativistic irreversible thermodynamics has had a rather pecuhar history. The 
first two proposals on the subject came from C. Eckart [l| in 1940 and from L. 
Landau and I. Lifshitz in the early fifties ^\. Since both of them have been shown to 
be particular cases of the so-called 'first order theories' ^, we concentrate here on the 
former proposal, mainly based on the construction of an energy stress tensor where 
heat flux is included, thus incorporating heat in the same status as a mechanical 
energy. Moreover it is also built so to satisfy the first law of thermodynamics. It 
was later shown by W. Israel and coworkers that both theories may lead to 

transport equations of the parabolic type, therefore violating causality. Also, it has 
been recently suggested that heat, being a non-mechanical form of energy, should 
not be included in the stress-energy tensor to be coupled to Einstein's tensor in 
the context of general relativity [6|. Regarding causality, we will leave additional 
comments on this interpretation and appropriate status in statistical physics for a 
future publication. Here we want to concentrate on the role of heat in deriving 
transport equations. 

Indeed, besides these difficulties, during the last two decades of the past cen- 
tury, several papers appeared pointing out another shortcoming of the two theories 
mentioned earlier namely, that they predict hydrodynamic instabilities in the linear 
regime [2, 13|, 15||. This interpretation is not correct. What they really show is that 
fluctuations under such conditions do not obey the Onsager's regression assumption. 
We remind the reader that this hypothesis states that spontaneous fluctuations that 
occur in a fluid in equilibrium evolve towards the equilibrium state following linear 
equations for the macroscopic variables [7[. This fact led workers in the held to con- 
sider the so-called 'second order theories' J, ^ as those which correctly solve the two 
main objections. This approach has now been taken as the basis of a new version 
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of a modified form for the relativistic hydrodynamic equations suitable to deal with 



the hot dense matter produced in the Relativistic Heavy Ion Collider, RHIC 

In this work we show that the problems mentioned above are caused by the struc- 
ture of the constitutive equation for the heat flux proposed by Eckart which couples 
it with the hydrodynamic acceleration. This relation, when introduced as a closure 
to the system of hydrodyamic equations, yields the so-called generic instabilities. It 
is important to point out that this coupling, whose thermodynamic meaning has 
always been questioned, has so far not been supported by experiment (lo| . 

As an alternative, we study the effect of a Fourier-like constitutive equation as 
can be obtained, for example, from Meixner's formulation of linear irreversible ther- 
modynamics jn], Q. This theory is based upon the idea that the heat flux is a 



13| and leads to a system of linearized 



component of a separate total energy flux 
transport equations which correctly describe the evolution of the fluctuations accord- 
ing to LIT. Moreover, the second law of thermodynamics is neatly derived and the 
stress-energy tensor maintains its canonical form including only state variables. 

These results suggest that those attempts seeking to deal with these problems bas- 
ing their ideas on Extended Irreversible Thermodynamics, EIT {l^, may be unnec- 
essary. The results obtained from A Meixner-like approach are free from additional 
parameters other than the standard transport coefficients. On the other hand, EIT, 
whose constitutive equations are of the Maxwell-Cattaneo type introduces additional 
parameters that, on the long run, have to be adjusted. Other difficulties with EIT 
theories have been extensively discussed in the literature [isl . 

To pursue our discussion we have structured the paper as follows. In section 
II the relativistic transport equations are derived from Eckart's proposal for the 
stress-energy tensor. Section III contains the linearized system and the consequent 
derivation of the so-called instability of the transverse velocity mode. Final remarks, 
including the alternative proposal corresponding to a Meixner-like scheme, are in- 



eluded in the last section of this work. 



II. RELATIVISTIC TRANSPORT EQUATIONS: ECKART'S FORMALISM 

The evolution of a relativistic fluid is described by the balance equations consid- 
ered together with suitable constitutive relations. The continuity equation in rela- 
tivistic hydrodynamics is a statement of conservation of particles. For a non-reactive 
fluid and denoting the particle flux as = nu^, such an equation reads 

Ni", = (1) 

where n is the particle number density. Here and throughout this work greek indices 
run from 1 to 4, latin indices from 1 to 3, and a semicolon indicates a covariant 
derivative. The fluid four-velocity is m^, so Eq. ([T]) can be written as 

h + nO = (2) 

where 6 = wf*, and a dot denotes a total time derivative. 

The general conservation equation for the stress-energy tensor 

Ti:.,, = (3) 

encompasses both energy and momentum balances. In Eckart's formalism, T'^'^ is 
given by 

ns 1 1 

= —u^'u, + phi;; + < + -g^u, + -u^'q, (4) 
c c c 

where e is the energy per particle in the comoving frame, p is the local pressure, tt^ 
is the Navier tensor, is the heat flux and /i^ is a spatial projector deflned, for a 
(+ + +—) signature, as 

K = sii + ^ (5) 

The last two terms in Eq. (Jll) are identifled as the heat flux contributions included 
in Eckart's formalism as part of the proposed relativistic generalization. Also, these 



terms satisfy orthogonality conditions equivalent to those for the viscous dissipation 
namely, 

uXu = = 0, = Q^^f^ = (6) 

It is worthwhile at this point to emphasize that these heat flux terms are not present 
in the stress tensor in the non-relativistic theory as well as in the relativistic version 

n 

of Meixner's thermodynamics |13l |. 

Equation (l3|), after computing the covariant derivative of T^^ given by Eq. ((41), 
using Eq. (I2|) and the operator ( ) yields the momentum balance equation 



+ ^ + ^''^-iM + ^1- + = (7) 

The evolution of the internal energy is given by the projection of Eq. ([3]): 

u'T!:.^, = (8) 
from which one can obtain the equation 

ne+pe + + q^^ + ^ii^q, = (9) 

where use has been made of the fact that, from Eq. (I6|), u^Qv^^ = —u'^f^qu and a similar 
relation holds for the viscous term. It is convenient to recast this equation in terms 
of the temperature as state variable instead of e by means of the local equilibrium 
assumption. That is, since e = e {n, T), one can write 

Using the relations (||)^ = + ^ and (y^) = C'n where 13 is the volume 

expansion coefficient, kt the isothermal compressibility and C„ the specific heat at 



constant n, Eq. ([9]) can be written as 

nCnf + (^^^ e + + gf^ + i^^g, = (11) 

Equations ([7]) and (fTB . the conservation equations, form an incomplete set. 
In order to obtain the dynamics of the state variables one has to introduce consti- 
tutive relations as closure equations. In order to propose these phenomenological 
expressions for the fluxes Tt^^ and , one calculates the entropy production from 
the entropy balance equation which is obtained by means of the local equilibrium 
assumption, that is, if the entropy per particle is s = s (n, e). 

Substitution of Eqs. ([2]) and (l9|) in Eq. (fT2ll and using the thermostatic relations 

(s^)^""^' (9^)^"r' ^^^^ 

yields an entropy density balance equation which can be brought to the following 
general structure 

ns + Jf,]., = a (14) 

where J^^-^ is identified as an entropy density flux and the entropy production is given 
by the expression 



a=--^— + -u,)- ^< (15) 

In order to assure the positiveness of cr, and thus satisfy the second law of thermo- 
dynamics, Eckart proposes the following constitutive relations 

< = -2r^Kh^,T^ - CeSi: (16) 

g^ = -Kh; (t'^ + Jtt^) (17) 
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where the transport coefficients involved are the bulk viscosity C, the shear viscosity rj 
and the thermal conductivity k. In Eq. (fT6|l r^^ is the traceless symmetric part of the 
velocity gradient tensor. Equation f flTl ) deserves a closer look. In it, the first term in 
parenthesis corresponds to the usual Fourier-type constitutive equation. The second 
term, which arises from the inclusion of the heat terms in the stress tensor, is not 
in the canonical form namely, it cannot be considered a thermodynamic force. This 
objection is equally applicable to the second term of Eq. (fTSll . A complete discussion 
of this issue can be found in Ref. ^ and will not be repeated here. However this fact 
is brought to the attention of the reader since, as will be shown in the next section, 
this term leads to the so-called instability found by Hiscock and Lindblom [^. 



III. LINEARIZED EQUATIONS FOR THE FLUCTUATIONS 

In this section, we shall study the behavior of the linearized equations which result 
from Eckart's scheme when spontaneous fluctuations of the state variables occur 
around equilibrium. These equations arise upon substitution of the constitutive 
equations given by Eqs. (11611171) into the general conservation equations (l2|), dT]) 
and f fTTl) . Next we linearize by setting T = Tq + ST, = 5u^ {5u^ = 0) and 
9 = 69, which according to Eq. ([2]) is equivalent to the condition n = uq + 5n. 
Here the naught subscripts characterize their equilibrium values and 6 denotes the 
corresponding fluctuation. The ensuing process requires a minimum effort. In fact, 
by simple inspection, one notices that the second bracket and the flrst three terms in 
the third bracket in Eq. ([7]) contain at least quadratic terms in the fluctuations as well 
as the third and flfth terms in Eq. f fTTll . Therefore, the only terms that introduce a 
difference with the results obtained in ordinary relativistic hydrodynamics are u'^Qu-a 
in the momentum balance and in the energy balance equations. Both arise from 
the inclusion of heat in the stress tensor. 
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The rest of the procedure is straightforward. Using the well known techniques 
of standard linearized non-relativistic hydrodynamics H, [3| we first write the 
linearized equations, 

5h + n^Se = (18) 

1 11 

— (no^o + Po) Siiu H (5n ^ + - — ST^ 

- 2rj (drr,) ^ - ^6T, - '^Su, =0 (19) 

nCjf + 66 -K [sT"" + J^m'-) = (20) 

where use has been made of the fact that p = p{n, T). 

Equations f|T8ll-fl20ll are the linearized equations for thermodynamic fluctuations 
in Eckart's version of relativistic hydrodynamics. The next step is to separate 69 from 
the transverse velocity. For this purpose we calculate the curl and the divergence of 
Eq. (fT9ll . The second operation yields 

rri -| 11 

- + ^ (^0^0 + Po) Se + —V^6n + -—VHT 

C + ^v) ^^Se - ^VHt = (21) 

This is a scalar differential equation with unknowns 60 (or 6n) and 6T . We shall 
come back to it afterwards. On the other hand, the first operation, recalling that 
the curl of gradient vanishes, yields 

kT 1 

-^6U, - - {noEo + Po) 6U, + 2r^WHU, = (22) 

where 6Ua = ^au^^^,^l ^"^^ the components of the curl of the velocity field, e^^ is the 
well known Levi-Civita's tensor. Thus, this procedure decouples the transverse mode 
fluctuations from the system of hydrodynamic equations, yielding an independent 
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equation. Calculating a Fourier-Laplace transform with transform variables and 
s respectively, we obtain 



5U4k', 0) ^s-^{noeo+Po) 



+ 



at 



t=0 



■s^ - 7^ (noEo +Po)s- 2r]k'^ 



(23) 



where dU^, is the Laplace- Fourier transform of 6U^. Thus, the exponential growth, 
or decay, of the transverse component of the velocitiy fluctuations will depend on 
the roots of the denominator, namely, the solution of the dispersion relation 

-To^._fnoeopo\^_^^j^,^^ (24) 



C4 V C2 



Equation ( I24I 1 has two real roots, namely 

' ' (25) 



C2 



2kTo 

which is precisely the result the authors of Ref. [3| arrived at in a more laborious 
way. Notice that, not only the procedure here developed yields the dispersion relation 
in a more clear, direct way but it also highlights two important points. Firstly, the 
equation for the transverse mode is completely decoupled from the system. Thus, it is 
unnecessary to make any assumptions about the direction of the velocity fluctuations 
as done in that work. The second and most important point we wish to emphasize 
is the fact that, as a consequence of decoupling of both components of the velocity 
held, longitudinal and transverse, one can easily identify the source of the exponential 
growth found in the transverse mode fluctuations. Indeed, for spatially homogeneous 
perturbations, k = 0, Eq. ( l25i l yields a positive root 

. = !^2^c^ (26) 

Equation (l26l) is a critical result. It clearly indicates that Onsager's assumption 
of the regression of fluctuations is violated. This is in open contradiction with the 
tenets of LIT. In this sense, Eckarts's formalism is suspect. 
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IV. DISCUSSION 



The results obtained in the previous section, have motivated the formulation and 
use of generalized constitutive equations. However, as can be clearly seen by inspec- 
tion of Eqs. (I21II24I) its cause is precisely the presence of the acceleration term in the 
constitutive equation for the heat flux, Eq. (fTTll which, in turn can be traced down to 
the phenomenological formalism flrst introduced by Eckart. The coupling proposed 



in Eq. (llTll is thus responsible for the theory being cathegorized as unphysical, and 
thus displaced by others jl^. 

On the other hand, an alternative proposal for the stress energy tensor, the one 
consistent with Meixner's theory can be easily shown to yield a system of hydrody- 
namic equations in which fluctuations behave canonically this instability is absent. 
Indeed, as discussed in Ref. [l3|, the stress energy tensor there assumed is 

T;: = pu^u, + ph; + n; (27) 

Equation (l27ll provides an alternative manner to introduce heat without identifying 
it as a state variable, consistently with classical thermodynamics. Heat is a form of 
energy, energy in transit, and must be included in a total energy balance equation. 
In Meixner's irreversible thermodynamics, heat flux is firstly introduced in a total 
energy balance given by 

^[E];. = (28) 

where Jj'^] is the total energy flux and includes the mechanic and internal energy 
fluxes as well as the dissipative heat flux 

Jj'^j = pc^u^" + neu" + q" . (29) 

This procedure yields an entropy production in terms of products of thermodynamic 
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forces and fluxes exclusively, 



T 

^ = -^^^ - fK (30) 



which is consistent with Clausius' idea of uncompensated heat and motivates a con- 
stitutive equation for the heat flux where the acceleration term is absent. Indeed, the 
constitutive relations in this formalism correspond to laws of proportionality between 
forces and fluxes which ultimately results in a momentum balance equation where 
the second derivative of the velocity is absent, and thus the evolution equation for 
the fluctuation of the velocity variable 5Uu reads 



C2 



5U^ = 2r]V'6U^ (31) 



which, performing a similar analysis as the one that lead to Eq. ( 1231 ) yields the 
dispersion relation 

^no.o+Po\^_2^p^Q (32) 



c 



,2 



Equation ( 1321) clearly predicts only an exponential decay in time for 6Uiy This result 
is compared with the unphysical behavior found in the previous section by analyzing 
the non-relativistic limit of both results for A; 7^ in Appendix A. As shown in this 
appendix, the root which generates the exponential growth mentioned above contains 
the thermal conductivity even in the non-relativistic limit. This is completely at odds 
with classical hydrodynamics and with the fact that the heat terms in Eckart's tensor, 
Eq. (JID, are assumed to be strictly relativistic. 

To finish this discussion, we would like to go back to Eqs. (1201 ) and ( I2T1) . a set 
of coupled equations for 6T and 69 (or 6n). To solve this system one has to go 
to the frequency and wave number representation and transform them into a set of 
two coupled algebraic equations. With the solutions obtained one can calculate the 
density-density or ternperature-temperature correlation functions. The former one, 



as well known 



3, 



17 



ISj, is related to the so-called dynamic structure factor whose 
11 



form, for a simple fluid, is known as the Rayleigh-Brillouin spectrum. It consists of a 
central, or Rayleigh peak and two symmetrically located peaks known as the Brillouin 
peaks. Rayleigh's peak has a width in frequency Acu cx Dfhk'^ where Dfh = is 
the thermal diffusivity and k = 27tX~^ when calculated with Eq. f l27ll . On the other 
hand, if Eckart's tensor is used, it appears the the Rayleigh-Brillouin expectrum does 
not exist. This fact would be in open contradiction with experiment. Calculations 
along this line will be published elsewhere. 



Appendix A: NON-RELATIVISTIC LIMIT 

In this appendix, the non-relativistic limit of Eqs. (l24l) and (l32l) are explored. In 
such limit, one can identify uqEo — > poc^ and further notice that po <C poc^ where 
Po is the equilibrium mass density. Introducing these facts in the dispersion relation 
obtained within Eckart's formalism, Eq. (I24ll . yields 

- pos - 2r,k' = (Al) 

while Eq. (l32l) yields the linear equation 

- pos - 2r]k'^ = (A2) 

Equation (lAll) has two roots which, using a binomial expansion for SuTo'qk'^ <C p^c"^, 
are 

si ~ '— (A3) 

S2-^ + -^ A4 
kTo Po 

Notice that the first root, Si, corresponds to the decaying behavior found within 
Meixner's approach, Eq. (IA2I 1. This behavior, in which the velocity fluctuations 
decay due to viscous effects is more physical than the one given by the second root 
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S2, which corresponds to the so-called instability referred to in Ref. [3|. This behavior 
of the velocity fluctuations is due to both viscous and thermal effects, a result which 
we insist, is untenable in this limit. 
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